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ABSTRACT 


A theox«tical sttiufy is tiade of th« Kicdhhoff th«o*y wcwstlc 
dumping in long cyltnderlcttl tv&es with tJw oh|ectlve of inclndtog th« 
corzoctions due to the velocity slip, the tes^rature 
teifHsmtuxe fluctuations. Explicit results have heen obtained for the 
attenuaticm coefficient and coapared with the »ec«it unpwhlia^®^ 
of Rott and Oberai, The «i*ily assuspticm sade is that the Stores layer 
thickness is wach smaller than the wave length A and this is os®d latter 
in the solutioa of the tiisadsnitie equations. 
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LXST OP Sl^BOLR 

ccitstants first used in Equation (6) 
sound velocity, (yRT^)^^^ 
specific heat at constant pressure 

constant to be eiq^ressed in tem of the nonentun accoanodation 
coefficient 

constant to be expressed in tem of the mmvgf aeconnodstion 
coefficient 

complex qismtitity first used in equation (d) 
imginaxy quantity, 

Bessel's fteetlon of the sero order 

Wsll tesfteratare fluctuatiim paraaeter, defined in page 7 

conplex propa|piti«t constmt 

attcoiuatien coefficient 

instantaneous acotattic pressure 

6as constmt 

inner radius of the tdbe 

instantaneous acoustic tmfMsratura 

instantaneous axial floe velocity 

iastmtaneous radial flow velocity 

Prandtl nusiber, iCpA 

Stohes aeem^ viscosity coefficient 

^mdc vistemity 

kinesatie viscosity 

dmsity 

fraquensy of osciHatim of the piston 
rstio of speeiAc hosts, Cp/C^ 



mm txm piith of the aolecules 
eonstentSf 4efln«d in page 6 
wave-length, m/m , 



I. INTRODUCTION 


The study of siaall acoustic disturbances into a nediuia in which the 
viscous and theraal dissipations are not negligible is essentially a Gasdynandcs 
Problon, Special interest attaches to this prcblea If the velocity slip end the 
ten^rature jump are not ignorable. This case is the subject natter of the 
present work. 

The classical theory of acoustic dashing in a long cylindcrical tube had 
its origin over a century ago with the work of Klrchhoff (1868), The above theory 
that takes into account viscosity at the tube wall is due to iiolniholtt (1363). The 
discrepancy of his theoretical results with the experlnents of Kundt (1968) Etotlvato 
Kirchlioff to present the conplete theory that includes the previously neglected 
effects of heat conduction in the fluid. Starting froa the linearised continuity, 
nonentua and energy equations, Kirchhoff obtained an algebraic equation for tho 
propagation constant in. He did this by assundng tho solution for all tho variables 
to be of exponential fora as regards tho axial coordinate and by satisfying tho 
boundary conditions of zero particle velocity and acoustic teatperatuxe at the tube 
wall. The solution of tho above algebraic equation gives both tho real and tho 
inaginary parts for the propagation constant is. While the real part accounts for 
tho attenuation factor for the waves progressing into the medium, the velocity of 
propagation of the disturbance can be derived from the imaginary port. 

The controversy arising from the disagroement of many experimental results 
with the classical theory is well known. Decades of scientific research have gone 
into tho study of this problem. The majority of tho researchers state that the 
actual effect of the tube is considerably greater than that given by Klrchho££*s 
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foYMula. Starting fro® the Kirc^off solution for th« acoustic quantities U, V 
«i(i T, it is possible to get few cxtx« teras for the attenuation coefficient. 

This is the starting point of the present investigation in which tiie corrections 
due to the velocity slip, the te^iexature j|uaf» and the tesfMirature fluctuations 
at the tube wall are Incorporated in the Krichhoff equations, 

Aamg the nuaerous invest lgati<ms directed towards finding the causes 
for the disagreewent of the Klrdihoff's results with the experinents, the work 
of Henry (19S1) is notable, He was the first to include the velocity slip, the 
teMpexature Juap and the toaperature fluctuaticms at the tidse wall in the 
Xirchhoff equations. His calculations were based on physical and aathenatical 
reasonings. He was aainly interested in the study of the ratio of specific heats 
of gsses by the sound velocity aethod. He incorporated various correctiwis in 
his analysis sepantely «mo by one and gave explicit foraulae for the attenuation 
coefficient and tho velocity of propagation. By neglecting tenw of the order 
JL in his calculations, he concludes that the velocity slip hss no effect in the 

attenuation coefficient but addb a snail tern to tho velocity of propagation. In 

tl 

addition, in tho detemlnant eaq>ansian, he nakes the approxinatiott a •“ in the 
terns containing the coefficient Such approxlnation effects the terw of the 

Qwmfw • 

The above iiwestigatlon wsts further extended by Weston (195S). tte i^wv 
the nathenatlcal spproxinatlons involved in the classical wcntk and included note 
terns in the Bessel functloii esqpanslons. The resistance - concept found in his 
work for the narww tubes is worth mmtloning. This concept energes from the 
uork of Lai* (l«9t). They have Included this concept in the analysis ftor the cas. 
when the tdbo radius is nuch snaller than the boundary l«y«t thldoiess. m also 
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••Iwt th# i^^ioxlitteion « « ^In hi# mslysis [cf. his wtimtiosi (56)|. 

tti« pi^per by Siields «t «1. (1065) pwsents th« xvsults of « mmorioil 
soliitioit of the Klrdl^ff oqusticois for tho piraf»sistiofi oqmstwnt «• thoy 
voport sociitste titawricsl vostilts for the sttemtstlen coofficlent end the 
velocity of piN^agetim derivable free the constant x* The roots Ij xaid 1^ 
are directly i^tsined by solving niiaesically the quadntic equation in 1. the 
first approjtlMite value of m (i.e,, <» ■ ^5 is incorpovifted in the argMent of 
the Bessel Ihnetioiis. Them the equation, obtained frwa the detendlnuit 
eapansion, is solved for a new value of n and is apdn substituted back into tlw 
equation as well as in the arfunent of the Bessel functions. The iteration ia 
continued until the required eecuracy is obtabied. they conpare the ntaaerical 
resttlta with tone eiqperisMnital aeasttreaents. Miile the attenuation factor shows 
good agraewmt, the velocity of propagatimi neods corrections duo to tho velocity 
slip and tho tn|>orature jvnp« particularly at low pressures. 

In stt unpubliidiod work, ftott and Obomi obtainod new aolutiona In series 
fom for the Kirdiboff theory of ecoustic deeping, they eaployed an espansion 
adhene with pavaneter a, wiiidt SMy be interpreted either es the ratio of the 
Stfdtea»lnyer thideneaa to tho wave Imgth or to the square root of tite ratio of 
the frequency of oscillatiena to that of collisions of the gas amleeules, to 
solve ttm psoblen and inoorpersted corrections due to tho velocity slip, the 
teiqierature Jnap and the wall teig>erature auetuations. Pvmi the goveming 
equations they have ebteined differential equations for the ecoustic pressure, 
the coaplete solutions of whidh give the desired results far the ettonuation 
ceefficiaat. the mto^ » the first • end the second « order terns are calculated 
separately, the anjor iinplificatiim of their erpsnsion schene is that tdiila 
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th* velocity slip and the teaperature jump affect the fi»t-order tors®, the 
tertt* owing to the radial pressure gradient and the axial dissipations do not 
enter till the second-oider of the achene. In the present worh. we shm# that if 
the velocity slip, the tesperature juap and the wall tcBperaturo fliM:ttsations 
be directly included In the equaticms originally fonulated by Xirdihoff, tiie 
results are identical to those of Rott and c»>erai. 



IX. OROaiOFP EQUATIONS AHD BOUNDARY CONDITIffllS 


We axe ecNticemeil with the acoustic oscillation of a perfect gas 
eontaiiied in a sesd-infiRite cylinderical tube, the gas is driven hy the 
periodic Motion of a piston located at the origin of the co-ordinate systen 
being enplc^ed. the ain of the present investigaticm is to obtain the 
attenuation coefficient of this snail disturbance. In the actual case, <me 
needs in the governing equations both the dissipative and the non-linear 
convective terns. Since the undistuid>ed gas is at rest, the velocity 
conponents are directly given by their acoustic quantities which tan considtexed 
to be infinitesinally snail] consequently, all the non-linear teras drop out. 

In the i^linderieal co-ordinate systen (x,r,6), the axl-syniwtric tine- 
dependwnt equation of continuity, the Kavier-Stcfltes equations for the nonentun 
and energy are t 
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( 2 ) 


( 3 ) 

( 4 ) 


Associated with the above four equations is the equation of state of the perfect 
(thomally as well as calorically) gas * 
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p • R p T (S) 

l€»f course, the mamaptim of « perfect ges is not essentisi} 

Polionring Ri:n:}thoi^, the e9q>xos8ion5 for the veloctty- coapoimits u 
an4 V mid the tenperature T, which ere essuned to iH^eve like a ftsnetl<m of r 


iMltlpliod by e"***^* , are given as 
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I 1 

T » (Y-l) (Aj . Xj,)^3 ♦ Aj [rCn^ - xp^]} 

(8) 


where »id X^ am the smll end large roots of the equation 



KlrcMioff ebtaiiMid the idwre solution using the cmly boundary cfuidition that 
along the axis of the ttdie, tlM various quantities reaain finite, the detailed 
derivation of these eqnatiens is glvwi in the hook theory of Sound Vol. II hf 
Lord Rsyloi#. 
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Boundary conditions : 

Th® propagation eonstmt a was detexsdned by Kirddioff by tising th« 
boimdary ccmditions of »«ro particle velocity and aoiustic teaperature at the 
tub® wall, Itowever, as stated in the tntroductlwi, ^® aotivation of the 
present woric is to incorporate the corrections due to the velocity slip, the 
temperature jinp and the tm^Mirature fluctuations at the tt&e wall in the 
Kirehhoff equations, the wall temperature fluctuations are governed by a 
paraaeter R, defined by ^ ahidi beciraes significant at low 

- VsS 

temparaturo d«« to T" depandenee of C, and ar® respectively the 

themal conductivity, the density and the speelHc heat of the tube mterial). 
Idiile, the phefKHMHia of slip and teiperature Jtrap are the part of rarefied 
gasdynaaics corresponding to slight dkparture froa the ecmtintaia Halt, these 
effects which ap{M«r as small corrections to the ccatiauua no slip, no 
tei^iierature |iap boundary cwiditions are significant when th# pressure and 
consequently density of tho gas Is low, th# corresponding eiqnressians for 
these corrections are given below t 

1} Wall tenperatture fluctuations 

We need to solve tho heat conthietimt equatiim in tho tidwi wall, Let 
j ^nx4>ht ^ uj |0 fiuetuating temperature at a point In th# tube aaterlal at a 
distance y fr«i th# imter surface of th# tube wall, Sim* Wig. 1. Two assu^ptiims 
axe Bade Here, naaely tho ojtlal heat ciswhictlon Is negligibl# in em^axiam to 
the radial heat conihietion and the inner surface is considered as plan# instead 
of cylindericali the snail penetration of the heat flux free tha gas into the 
wall provides justification tor these asauivtions (cf. Ref, 2) 


We llien have 
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with th« boundaTy C0B(B.tiims • 0 , y - and « T|^ , y « 0 
whew r^f 1$ the instantaneous wall teaperatuw. 
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The solution of the equation (iO) ®ay be obtained as 

/$ ^s **^1/2 

-C-V=— ) ' y 


Tg - e 


( 12 ) 


Making a heat balance along the inner surface of the ttd>e wall* we have 


vHl . k !!il 


(IS) 


After sinplifieatim, we obtain 


« . K (|.)2 H 


( 14 ) 




ii) Tenqperature juiqt 


The statement of the teafierattire jtaqp at the tube wall is as follows 
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On substituting fbr fro* the equation (14), the following boundary condition 


is obtained t 
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ill) Slip 

TImi pmseribed bouaii«x]r conditions «t tho surfsee of iwill «r© 
tb?t tho noxMil velocity is toxo; hoifevor, dtte to slip tho tengentiel velocity 
is pxv^Ttionsl to tie gradient 0 . Explicitly, 
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III. METTiOD OP SOUmCN 


the qmJlltfttlvtt m well as qumtltative affects ©f the factots 
4esarlbe<! In the pvevlims chapter can be i^tained by using the notified 
boundary conditions derived In that ehqiter. The corrections are sought in 
the various acoustic t|iiatitities u, v and T. The corresponding expressions 
after satisiying the boundary conditions at the tub® vail will give 
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( 21 ) 


Kriting the# In ^ natrix fom with A^, Aj, and as colian vectors, we 
f*t PI CAl • C01« 

coeffi«i«Bt aaitrix M contains the tnknowns X|, Xj imd tlwi paranater a. 
the aethod of obtaining the propagation constant n is siaHar to 
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that of Iwirdihoff. lha emtdltion that the deterwlnant of the coefficient 
nattix should vanish for the ncm>sero valws of the vector [Aj, A 2 » Aj] gives 
an algohraic equaticm for the propagatlem ccsistant. In interpreting the 
solution there are two special cases of Interest, the first, id.de tidje cwte 
is the saiae as that considered by KlrdihoK. The e 3 q»lanaticHis given by 
Kirdtiioff showed that the approxiaate values of X| and ■ are both equal to 

2 1/2 
!L- . froa dlls he concludes that the tew *rC*^ - X,) ' * appearing as an 
a2 

arguaent in one of the Bessel functions is to be reganted as sMilI. However, 
the other two ar^saents 'rCa^ - end *r(«^ - X 2 )^^^* Mte large for diis 

ease. In the other oxtreae case of narrow tube, all the three arguewits are 
to be regarded as saall. The relevant asyaptotic expansions of the Bessel 
functien are given below s 

a) IWien |lj is large and such that its ieaginary part is positive 


|yfn • -1 Cl ♦ 


-.1 • U 


When iz| is smll 
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Having diio»s#d the aj^roxiwtioas involved In the series 
expansions of the sere order Bessel function, the next task Is to obtain the 
roots Xj and Xj. These roots, appearing in the Klrdrfioff oquatlons, are 
independent of die tdbe radius, being detewlned by the equation 


BX^-CX*»»0 


(25) 
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a»iving th« qita4initle eqimtion (2S), we get 


C26) 


^1*2 21iL* hi^er owlet teww] 


(27) 


III the hinoniel exi^ension for the expression under the square root sifa, we 
have used the fact, « i. (prom the fona of the expressions for B, C and 
0, one »ay show that this condition is i^ilfilled if the Stidies layer thickness 
is Mich Simller than the wave length). 

Sidistitutlng for 1, C, and D in the equatiem (27) and after 8ls^>lifi cation, 
we obtain 
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i) hlDi! TUBE 

In the previous section we have jaeiitioned the ICirchtuiff i^iproxiBiations 
in the Bessel function expansion with regard to its arniMmt. Our ale now would 
be to include these iqq>rosd«atio»s in the Kirdihoff detersdnent (diter»ittant 
of the coefficient aatrl* M) before eiqianding it. Bor the present case, we 
need the asyaptotic eiqiansions of the sere order Bessel ftmctlon given by 
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equations (23) an^ (24). in the es^wslons given belov «o have xetained two 
terws Instead of one as given in Rayl®i^»s hook. A ctMi^ariscm of the results 
of the present Investigation with the recent unpuhliahed work of Rott imd 
Oboral is not possible without the second term appearing in all the throe 
exnimsioas. ilence we have 




(30) 


Incorporating these erpretsiont in the coefficient natrix M, we expuid its 
deteminwit. Ihls* alter divisium by Q gives 
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By ttsiaf the approxlwatloiw (2S) « (29) for Xj « Xj «id siapliiyiaf we obtain 
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Th« next talk mnli be to obtain the quartic equation end tliwa. Its solution* 
These ave as foilcwis i 
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Out of th« tiio |KMisil»l« toot* glv«fi by oqiiwtiim (33), we bwre eJiosett the tsa® 
whic^ beconet Identical with the value of obtalnedi by neglectinf viscosity 
and theraal conductivity la th® foveming ixiiwtioits. Putting h • iw (so that 
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1 ? 1 

(U*) « (|j-) (l-i) in t^e equation (34), v® get 
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th. I,.i ,,.rt ft,, ptop,B.tlon constmt mi 

0 . iiig tlio ne^ivc sign (comsponds to das^iing), wo have tlio following 

c 3 c:^iKss.^iojj for the attenuation ooefiicimt 
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ii) mmm tube 

Ttiis case arises iditn the layer affected by friction, instead of foisting 

a thin coating to the valla, extends itself over the tdiole section. This 

1/2 

hofpens when the f!tokes«liQwr thickness, ^fined by (^> ^ , is eimarahle to 

the tube radius. In view of th# assixaittion * the arginsent of all 

Bessel fBnctiona in the Matrix H, becoraos snail. Hence, a fresh start is 
necessary for the yieswit aime. It way be pointed out that th# aathod diseaissed 
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by Rott md Ob«r«i u»4m to toaalti valid for arbitrary ratios of Stokas- 
layar thleknoss to tha tidM radius. In th« osqtrrsnlous givcm below, we 
have to retain three ton» of the ejiir^isioits In order to coajpare present 
results with those of Rott and Clierai. fJonce, we have 


& 

ll« In Q « 
In Q| • 

Q 2 * 


f C»^ - ^ II » ^ t«^ - ^ ^ 1 


t /_2 


J •» 


(wT . k^m ♦ f- *^iwir • k{)^ ♦ .. 3 


I (»^ - xptl ♦ ^ in - V ♦ ^ C»^ - ♦ •• 3 


(37) 


After Ineorporating these expressions in the eoefflelent natrlx M, we e^qpaad 

iu Tht., •tut .lupllflctl™ mi division by Q Qj Qj glvos 


mZ md 


(* 




48 


o .2 


IC«" 


r^ r^ 

WmM m 


• Ij) ♦ -? {•* - ^ ** 


Fv r^ I, ^ ^ ^\4 % 


2a 


4a* 


rEh rf 


*w *11 4 2 t »#. 

{1 ♦ « «2 - rrr-» rl tU 


2 .2 

*iyf M 


I ; ^ t» Is2 ^2 |,.2 

t * 11 ^ ^ ^ ^ 


fw 48 


24 V 


48 


2a 




4a* 


0 a"^ a^ 


F^v^rix^ 

4 a^ 


>3 • 0 C3«) 


f)tt ftirther siapliflcetlon, we obtain 
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,]4'' 

^24 




II h 


2m 


)] • ‘ l—‘-^ 


»; Hyl) vh 

“SS 5 ? 



* h 

a 



(1 ♦ « 



Cl 


Pv tk >J 4 t| 


>1-0 


(503 


SolviniK till* qusrtic aqiiation aii4 dioosiRg t!i© pfoper loot, we get 


*2 . 0 vh 

M • 




4E w ^ IftE 

rT a 


2 Y 


2^2 P ohCt- 1) rj h^CY-l) 


W 

4.2 


«. (w ♦ «)} 

ar 


(40) 


by using the appsexiisations (28) I (29) toy & 12 * Sinplifying, we obtain 


1 

2C2wIiy)^ ^ 2P ^ ^2 V ©(y- Dh P^ tJ o^ h*(Y-13 

» • ♦ I , . . CC • ♦ ' “"j-.. *-y “ 

** Yw a % * ti ^ 8a^ y 

w w 




(41) 


Putting h • to in the above e«taation, we get 

1 


u i* ♦ 


2(l*i) (V«r)* asv 6E^v^ Fr,<.{T-l)l» 

(1 • — ♦ ' — 


a % 


% a a^ 


4a 


8 a-' 


^ C|a«)3 


im 



Ftoh thlSf «*pT«s*ion for tho ottonuation coofficlont^ corresponds 

to dawning. Is obtained 
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E^Ct)^ + |sl 


FoCy-1) 

2 



4(y)1/2 



1 _ 

(|^ «) y(y)^1 


i 

where |n^l • 


C4S) 



IV DISgtSSIOH 


m shall suMviss th® msults t^tained la th« prwious dieter 

•ad gl» e<»|»ri.„, wh«.«.r poMlkta. .1th tho» obtain..! b, .»h.t .nhoro. 

i) WIMJ TOE 

the sfftct of thtt ttall teaponttitre fluettaitloa paraaetor K is to 
docmaso tha thonal coBtrtbutloa to das,>ing. TO slip has a aogatlvo 
contrihutlost to dai^lng. 7h# sl|pi of tha eoatxibotloa of the taapaxatuxa 
dapands on mgnituda of K valatlva to 1). Soo Tahla 1. liansy and 
lifaston hami studlad tha offset of slip and toafteratuze jusp ai^ concluded that 
they have no affact on danplng. In fact, these authors did not cany out the 
Bossal ftmetion aspansloaa to the ozder i*- (as also pointed out by Rott and 

Tm 

(Utaxai)} this has bean dona in tha present tmxic. Ikmairwr, the next higher 
order tans witli the ooafficiant ja^j agree aith those of boston. In Pig. 2 
we amp&ta our results with the ejqiNsrlaeatal ■easuxenents reported in ref . 4 
for Argon. 

ii) nmm toe 

TO slip has a large and negative contribution to deifying. This 
result agrees with tdiat of Weston. TO <»nitributioit of the tos^rature juap is 
positive but of Mich snallor order. 

1hou|h> present nethod involves wre aathenatical calinilations than 
that of Henry end Weston, the generality thet no i^pproxittation of TO type assuMd 
hy these authors has been tiled Imtifies it. As stated earlier, presTO results 
are identical with TOm of Itott and Cborai. Finally, It nay be pointed out that 
TO contributieo of TO tern * ^ appeariiig in TO root alth 

g # 

TO coefficient (|H eases nad hmicse caimot be noglei^ed in TO analysis, 
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FIG,:2: SOUND ABSORPTION IN Ar 
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